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Abstract We perform a kinematic analysis of the Hipparcos and TRC proper motions 
of stars by using a linear Ogorodnikov-Milne model. All of the distant (r > 0.2 kpc) stars 
of the Hipparcos catalog have been found to rotate around the Galactic y axis with an 
angular velocity of = —0.36 ± 0.09 mas yr“^. One of the causes of this rotation may be 
an uncertainty in the lunisolar precession constant adopted when constructing the ICRS. In 
this case, the correction to the lAU (1976) lunisolar precession constant in longitude is shown 
to be Api = —3.26 ± 0.10 mas yr“^. Based on the TRC catalog, we have determined the 
mean Oort constants: A = 14.9 ±1.0 and B = —10.8 ±0.3 km s“^ kpc“^. The component of 
the model that describes the rotation of all TRC stars around the Galactic y axis is nonzero 
for all magnitudes, = —0.86 ± 0.11 mas yr“^. 
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INTRODUCTION 

Based on the Ogorodnikov-Milne model and using the proper motions of stars, Clube 
(1972,1973), du Mont (1977,1978), and Miyamoto and Soma (1993) modeled the Galactic ro¬ 
tation. This modeling showed that, apart from the Galactic rotation parameters, it is possible 
to control the inertiality of the catalogs being analyzed and to refine the adopted precession 
constant. Here, instead of the currently popular method for solving the Ogorodnikov-Milne 
model equations using parallaxes, we use a method for determining the kinematic parameters 
of the Galactic rotation that is completely free from parallax errors to analyze Hipparcos 
data (ESA 1997). In this method, the parallactic factor is assumed to be equal to unity. 
First, this assumption makes it possible to analyze stars even with negative parallaxes (the 
Hipparcos catalog). Second, since the proper motions of TRC stars (Hpg et al. 1998) can be 
analyzed in full only by this method, this approach yields comparable data for both ICRS 
catalogs (International Celestial Reference System). The kinematic method for controlling 
the inertiality of the catalogs of stellar proper motions is based on the analysis of the two 
components of the rigid-rotation tensor that describe the rotation about the Galactic y and 
X axes in the Galactic coordinate system. 

BASIC EQUATIONS 

In this paper, we use a rectangular Galactic coordinate system with the axes directed 
away from the observer toward the Galactic center (/ = 0°, 6 = 0°, the x axis), along the 
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Galactic rotation (/ = 90°, & = 0°, the y axis), and toward the North Galactic Pole (6 = 90°, 
the z-axis). In the Ogorodnikov-Milne model, we use the notation introduced by Glube 
(1972, 1973) and employed by du Mont (1977, 1978). When using only the stellar proper 
motions, one of the diagonal terms of the deformation matrix is known (Ogorodnikov 1965) 
to remain indeterminate. It is possible to determine the differences between the diagonal 
elements of the deformation matrix, for example, in the form {Mii — M^) and (M^ — M^). 
In this approach, the basic equations can be written as 

/i; cos 6 = (l/r)(X 0 sin/— Pq cos/) — ( 1 ) 

—Mj” cos I sin b — M" sin I sin b + M~^ cos b+ 

+Mj^ cos 2 / cos b — sin / sin b+ 

cos I sin b— 

—0.5(M+ — Mj^) sin2/ cosb, 

/ife = (l/r)(XQ cos / sin b + Yq sin Isinb — ( 2 ) 

—Zq cos b) + Mg" sin / — M" cos I— 

—0.5Mj^ sin 21 sin 2b + Mj^ cos / cos 2b+ 
sin I cos 2b— 

—0.5(Mj^ — Mj^) cos^ / sin 2b+ 

+ 0 . 5 (M 3+3 -M+)sin2/>, 

where Xq, Yq, and Zq are the velocity components of the peculiar solar motion, and M^, 
M^j, and Mg^ are the vector components of the rigid rotation of an infinitesimal solar neigh¬ 
borhood around the corresponding axes. In accordance with the adopted rectangular coor¬ 
dinate system, the following rotations are positive: from axis 1 to axis 2 , from axis 2 to axis 
3, and from axis 3 to axis 1. The quantity M^ is an analogue of the Oort constant B. Each 
of the quantities M{ 2 , M^g, and M^ describes the deformation in the corresponding plane. 
The quantity M {2 is an analogue of the Oort constant A. The diagonal components of the 
deformation tensor M{g, M^, and M^ describe the overall contraction or expansion of the 
entire stellar system. Equations (l)-(2) contain eleven unknowns that can be determined by 
the least-squares method. The quantity 1/r is the parallactic factor, which is assumed to be 
equal to unity. In this case, the stars are referred to a unit sphere. In this approach, all of 
the parameters being determined are proportional to the heliocentric distance of the stellar 
centroid under consideration and are expressed in the same units as the stellar proper motion 
components, i.e., in mas yr"b Before the appearance of the Hipparcos catalog, researchers 
were forced to use this method of analysis because of the lack of highly accurate stellar 
parallaxes. When the distances to stars are known, the parallactic factor is 1/r = 7 r/ 4 . 74 , 
where tt is the parallax, and the factor 4.74 is equal to the ratio of the number of kilometers 
in an astronomical unit to the number of seconds in a tropical year. To express the solar 
velocity components (mas yr“^) in km s"^, they must be multiplied by 4 . 74 / 7 r; to express 
any of the derived (in mas yr"^) components of the deformation and rotation tensors in km 
s"^ kpc"^, they must b e multiplied only by the proportionality factor 4.74. 

THE SYSTEM OE THE HIPPARCOS CATALOG 

We divided the Hipparcos stars into seven groups, depending on the heliocentric distance 
(kpc): 0.05 — 0.1, 0.1 —0.2, 0.2 —0.3, 0.3 — 0.45, 0.45 — 0.66, > 0.66, and the group of stars with 
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Table 1. Kiiiematio parameters inferred from the proper motions of Hipparcos stars 


Parameter 

1 

2 

3 

4 

5 

6 

7 

8 

r, kpc 

0.074 

0.142 

0.243 

0.361 

0.570 

1.2 

r ^2 

> 0.2 

N * 

13453 

29378 

20032 

16040 

11441 

10015 

3833 

58675 

A © 

25.32 ± 0.70 

14.27 ± 0.25 

8.10 ± 0.17 

6.05 ± 0.14 

4.53 ± 0.13 

3 .. 56 ± 0.11 

2.23 ± 0.19 

5.89 ± 0.08 

V © 

47.97 ± 0.69 

23.24 ± 0.13 

13.37 ± 0.18 

10.24 ± 0.15 

8.06 ± 0.14 

6.75 ± 0.13 

5.81 ± 0.18 

10.22 ± 0.08 

Z& 

18..59 ± 0.68 

10.04 ± 0.24 

5.66 ± 0.16 

4.27 ± 0.13 

3.23 ± 0.12 

2.49 ± 0.10 

2.03 ± 0.16 

4.10 ± 0.07 

' 0 

20.1 ± 0.2 

19.6 ± 0.2 

19.1 ± 0.2 

21.6 ± 0.2 

26.4 ± 0.4 

45.7 ± 0.7 

~ 62.0 ± 1.7 

22.0 ± 0.1 

r © 

62.2 ± 0.7 

. 58.6 ± 0.5 

. 58.8 ± 0.6 

. 59.4 ± 0.7 

60.6 ± 0.8 

62.2 ± 0.9 

69.0 ± 1.7 

60.0 ± 0.4 

B © 

18.9 ± 0.7 

20.2 ± 0.5 

19.9 ± 0.6 

19.8 ± 0.6 

19.3 ± 0.7 

18.0 ± 0.9 

18.1 ± 1.4 

19.2 ± 0.3 


2.40 ± 0.86 

2 . 20 ± 0.:30 

2.96 ± 0.21 

2.85 ± 0.17 

2.77 ± 0.16 

2.91 ± 0.14 

2.48 ± 0.21 

2.90 ± 0.09 


- 0.18 ± 0.70 

0.06 ± 0.25 

- 0.15 ± 0.18 

- 0.10 ± 0.15 

0.07 ± 0.14 

- 0.11 ± 0.13 

- 0.53 ± 0.21 

- 0.16 ± 0.08 


- 2.44 ± 0.69 

- 1.06 ± 0.25 

- 0.62 ± 0.18 

- 0.15 ± 0.15 

- 0.24 ± 0.15 

- 0.18 ± 0.13 

- 0.68 ± 0.21 

- 0.36 ± 0.09 


- 2.40 ± 0.68 

- 2.47 ± 0.24 

- 3.01 ± 0.16 

- 3.03 ± 0.13 

- 3.00 ± 0.12 

- 2.85 ± 0.10 

- 2.76 ± 0.16 

- 2.93 ± 0.07 

M +- M + 

0.08 ± 1.76 

- 2.13 ± 0.62 

- 2.00 ± 0.42 

- 1.40 ± 0.34 

- 1.16 ± 0.32 

- 1 .. 50 ± 0.27 

- 0.62 ± 0.41 

- 1.35 ± 0.19 


- 4.01 ± 0.92 

- 0..56 ± 0.33 

- 0.10 ± 0.23 

- 0.12 ± 0.20 

- 0.14 ± 0.19 

- 0.07 ± 0.17 

- 0.75 ± 0.25 

- 0.11 ± 0.11 

M + 

- 1.38 ± 0.92 

- 0..58 ± 0.33 

0.41 ± 0.23 

- 0.03 ± 0.19 

- 0.20 ± 0.18 

0.14 ± 0.16 

0.49 ± 0.25 

0.13 ± 0.10 


- 1 . 6.5 ± 1.73 

- 1.33 ± 0.63 

- 0.61 ± 0.44 

0 . 17 ± 0.:37 

- 0.07 ± 0.36 

- 0.24 ± 0.33 

- 0..39 ± 0..52 

- 0.14 ± 0.21 

Ixy 

-1 ± 11 

13 ±4 

9 ±2 

7 ± 2 

6±2 

7 ± 1 

4±2 

7±1 


Note: r is the mean distance, is the number of stars, the coordinates of the solar apex L© and B© and the vertex hy are in degrees, 
t'© is the solar velocity in km s“^; all of the remaining quantities are in mas yr“T 


negative parallaxes. We imposed a constraint on the stellar space velocity that allows ns to 
discard only those stars whose space velocities exceed the hyperbolic velocity, for example, 
stars with enormous peculiar velocities acquired through explosions and encounters, V = 
220 + 62 = 282 km s“^. Here, 220 km s“^ is the circular velocity of Galactic rotation at the 
solar Galactocentric distance (recommended by lAU 1986) and 62 km s“^ is the mean stellar 
space velocity dispersion in the solar neighborhood (estimated by Oort). The Hipparcos stars 
whose space velocities exceed the hyperbolic velocity were found, for example, by Moffat et 
al. (1998). In this case, when using only the stellar proper motions, we estimated the stellar 
space velocities on the basis of the Kleiber theorem (Agekyan et al. 1962): 

3 14 

|Ki|-4.74-r<^-|K|. 

Here, Vt = cosb + lA the tangential velocity of the star in mas yr“^, and V is the 
total space velocity of the star in km s“^. Table 1 presents the results of the simultaneous 
solution of Eqs. (1) and (2). Columns 1-7 give the parameters obtained from the proper 
motions of the stars located in seven “thin” spherical shells. All of the stars with negative 
parallaxes were assumed to be located at a heliocentric distance of 2 kpc. The last column 
gives the kinematic parameters determined from the proper motions of 58675 distant stars 
located farther than 0.2 kpc from the Sun (the mean heliocentric distance of this group of 
stars is 0.371 kpc). We emphasize that these are special solutions; they were obtained by 
setting the parallactic factor equal to 1/r = 1. We needed the parallaxes to divide the stars 
into distance-limited groups. This method made it possible, first, to use the most distant 
Hipparcos stars with negative parallaxes, second, to completely eliminate all of the effects of 


3 




Fig. I. Component of the rigid-rotation tensor for 
Hipparcos stars against stellar heliocentric distance r. 


random parallax errors when solving Eqs. (1) and (2), and, third, to ensnre that the resnlts 
obtained by this method are in close agreement with those obtained from the proper motions 
of TRC stars. Since we used this approach to control the inertiality of the Hipparcos catalog, 
the solution from the last column of Table 1 that is of greatest interest in this respect may 
be called inertial. Let us consider the derived kinematic parameters of the inertial solution 
(in mas yr“^): 

M+ = +2.90 ± 0.09, 

M- = -2.93 ± 0.07, 

M +3 = - 0.11 ± 0 . 11 , 

M- = -0.36 ± 0.09, (3) 

M +3 = +0.13 ± 0.10, 

M 3 - = -0.16 ± 0.08, 

M+-M+ = -1.35 + 0.19, 

M 3+3 - M+ = -0.14 ± 0.21. 

As we can see, the following quantities are statistically significant: M{ 2 , and (M{i—M^). 
Based on these quantities, we find that A = 13.7 ± 0.4 km s”^ kpc”^, B = —13.9 ± 0.3 km 
s”^ kpc-\ C = —3.2 ± 0.5 km s”^ kpc”^ and the vertex deviation = 7 ± 1°. To this end, 
we used the standard relations A = M 12 -4.74, B = MA -4.74, C = 4.74- (M{i — M^)/2 and 
tan24y = —CjA. 

The Oort constants A and B obtained by this method are in good agreement with the 
values recommended by the lAU (1986) and with the values from the paper by Bobylev 
(2004): A = 13.7 ± 0.6 km s”^ kpc-\ B = —12.9 ± 0.4 km s”^ kpc”^ in which they were 
calculated by analyzing their recent determinations by various authors. The vertex devia- 
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Fig. 2. Kinematic parameters inferred from the proper motions of TRC stars versus magnitude (Tycho B mag.). 


tion also agrees with its published determinations based on various observational data, for 
example, with those by Dehnen and Binney (1998). Note, nevertheless, that is also sta¬ 
tistically significant. Figure 1 shows a plot of against heliocentric distance constructed 
from the data of Table 1. As we can see from Table 1 and Fig. 1, is nonzero for the 
most distant Hipparcos stars. It should also be emphasized that, as can be seen from Fig. 1, 
the constraint r > 0.2 kpc used to obtain the inertial solution and to determine has 
such an effect that we obtain a lower limit on 

THE SYSTEM OE THE TRC CATALOG 

The results of our solution of Eqs. (1) and (2) based on TRC stars of mixed spectral 
composition are plotted against magnitude in Eig. 2. The random errors of all of the sought 
unknowns are small in all of the magnitude intervals, each of which contains ~ 100000 stars, 
and are equal to 0.1-0.2 mas yr“^ for faint stars. As we can see from Eig. 2, the parameters 
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Table 2. Precession eorreelif)ns Api and AE. [iias yr ^ 


Reference 

Data 

Api 

AE 

Miyamoto and Soma 11993) 

.4CRS 

— 2 . 7 ( 0 . 3 ) 

- 

Waller and Ma (1994) 

VLBI 

-3.6(i.i) 

- 

Chariot et ul. (1995) 

VLBI+LLR 

“3.0(0 2) 

- 

Rybka(1995) 

PPM 

“3.1(0.2) 

— 1-3(0.2) 

Bobylev (1997) 

PUL2-PPM 

— 2 . 8 ( 0 . 8 ) 

- 

Ma et ul. (1998) 

VLBI 

“2.84(0.04) 

- 

Vityazevt 1999) 

CGC-HIP 

“3.4(1.0) 

“3.3(1.0) 


that describe the deformation in the yz plane and the rotation around the x axis, i.e., M 23 
and are virtually equal to zero. The parameters that describe the deformations in the 
xy and yz planes and the rotation around the z and y axes differ significantly from zero. The 
values of are nonzero for the faintest TRC stars. We calculated the mean, 

M” = —0.86 ± 0.11 mas yr“^, (4) 

by using the results obtained in nine intervals of magnitudes fainter than 9"^.5. This result 
must be compared with the results of our analysis of the proper motions for Hipparcos stars. 
The deviation of the vertex Ixy is virtually equal to zero in each magnitude interval, which 
does not confirm the conclusions drawn from our analysis of the Hipparcos stellar proper 
motions. This result probably stems from the fact that to compile the reference catalog that 
the TRC was intended to be (Kuzmin et al. 1999), the stars were selected kinematically; no 
high proper-motion stars were used. We calculated the mean values of = 3.15 ± 0.20 
mas yr“^ and = —2.27 ± 0.06 mas yr“^ by the mean of the nine values obtained in 
the intervals of magnitudes fainter than 9”^.5 (the mean being m = 11™.2). This yields 
the following Oort constants: A = 14.93 ± 0.97 km s“^ kpc“^ and B = —10.77 ± 0.31 
km s“^ kpc“^, which are in good agreement with those calculated by Oiling and Dehnen 
(1999) from the Tycho/ACT proper motions of red giants: A = 14.2 ± 1 km s“^ kpc“^ and 
B = —12.7 ± 1 km s“^ kpc“^. The model component that describes the rotation about the 
Galactic y axis inferred here from a kinematic analysis of the proper motions of TRC stars 
is My = —0.86 ±0.11 mas yr“^ and confirms the rotation that we found from the Hipparcos 
catalog (solution (3)). In this case, we have an upper limit, since the data for the faintest 
stars include the effect of the real rotation of nearby stars that to the Local stellar system 
(Fig. 1). 

CORRECTION TO THE lAU (1976) PRECESSION CONSTANT 

The proper motions of the Hipparcos stars are free from the effect of precession of the 
Earth’s axis. However, the Hipparcos catalog is an extension to the optical range of the ICRE 
system (Ma et al. 1998), which is based on ground-based VLBI observations of radio sources. 
To determine the precession parameters from the analysis of My, we use two equations that, 
given the numerical values for epoch J2000.0 (Perryman et al. 1997), can be written as 

M^ = n^- 0.0965Api ± 0.4838AE, (5) 

My = ny + 0.8623Api - 0.7470AE. (6) 
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Here, f2a; and Q:y are the vector components of the real rigid rotation of the stellar system, 
and this rotation actually exists at heliocentric distances up to r ^ 0.2 kpc. Assuming that 
at large distances, Mj, — ffa, = 0 and My — Qy = —0.36 ± 0.09 mas yr“^ includes only the 
precession quantities, we obtain from the solution of Eqs. (5) and (6) 

Api = —0.50 ± 0.13 mas yr“^, 

AE = —0.10 ± 0.02 mas yr“h 

Setting AE = 0, we obtain the following estimate from Eq. (6): 

Api = —0.42 ± 0.10 mas yr“^. (9) 

Table 2 gives the results of the determinations of two precession parameters by various 
authors: the correction to the lAU (1976) constant of lunisolar precession in longitude, 
Api, and AE, which is the sum of the corrections to the rate of planetary precession and 
the motion of the zero point of right ascensions. Miyamoto and Soma (1994) determined 
the precession correction based on a kinematic analysis of the proper motions of ACRS 
stars using the Ogorodnikov-Milne kinematic model. Walter and Ma (1994) determined 
the precession correction based on VLBI observations of extragalactic radio sources from an 
annual catalog. Chariot et al. (1995) determined the precession correction by analyzing a 24- 
yr-long series of lunar laser radar observations and a 16-yr-long series of radiointerferometric 
observations. Rybka (1995) determined the precession parameters based on a kinematic 
analysis of the proper motions of PPM stars using the Ogorodnikov-Milne kinematic model. 
Bobylev (1997) determined the precession correction by analyzing the differences between 
the Pulkovo absolute proper motions and PPM. To obtain the ICRE system that is based 
on ground-based VLBI observations of radio sources, Ma et al. (1998) adopted the following 
correction to the lAU (1976) lunisolar precession constant: Api = —2.84 ± 0.04 mas yr“^. 
Vityazev (1999) determined the precession parameters by comparing the proper motions 
of CGC and Hipparcos stars. Our kinematic analysis of the proper motions of Hipparcos 
stars yielded Api = —0.42 ± 0.10 mas yr”^ (solution (9)). It may be assumed that Ma et 
al. (1998) slightly underestimated Api. Since the ICRS was constructed using precisely this 
value, we obtain a statistically significant “addition” to this correction, Api, when analyzing 
the Hipparcos and TRC catalogs. Therefore, the correction to the lAU (1976) lunisolar 
precession constant is 


( 7 ) 

( 8 ) 


Api =—3.26 ± 0.10 mas yr L (10) 

This value agrees well with the data given in Table 2. 

THE MOTION OE DISTANT OB STARS 

Results (7)-(9) have yet another important implication: they show that My = —0.36 ± 
0.09 mas yr“^ may have a nature unrelated to the real stellar motions. The H I layer in 
the Galaxy is known to have an extended warp (Kulikovskii 1985; Carroll and Ostly 1996). 
Smart et al. (1997) and Drimmel et al. (2000) considered the hypothesis that the extended 
Galactic H I layer affects the kinematics of OB stars. Drimmel et al. (2000) used simulations 
to estimate the angular velocity of the precession of OB stars in the zx plane, -25 km s“^ 
kpc“b The effect of our = —0.36 mas yr“^ (solution (3)) on the kinematic parameters 
of stars is easy to calculate. Since ■ 4.74 = —1.7 km s“^ kpc“^, it is equal to —1.7 and 
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Fig. 3. Linear velocity vectors = U' versus helio¬ 
centric distance r. The Galactic plane coincides with the 
xz plane, the x axis coincides with the r direction, and the 
z axis is directed upward. 



Fig. 4. Stellar residual velocity components It' versus 
heliocentric distance R. The solar motion with respect to 
the LSR with velocities {U. W IF)© = (9. .5. 7) krn s“' 
w'as taken into account. The dotted line indicates our 
dependence, which must be eliminated. 


—3.4 km at r = 1 and 2 kpc, respectively. The presence of snch a fictitions wave conld 
completely explain the precession velocity of -25 km s“^ kpc“^ fonnd by Drimmel et al. 
(2000) for distant OB stars. In Fig. 3, the linear rotation velocity vectors M?13r are plotted 
against heliocentric distance r. This rotation is attributable to the = —0.36 mas yr“^ 
that we inferred from the Hipparcos catalog. The directions of the vectors are indicated in 
accordance with the fact that the rotation takes place in the zx plane and that, according 
to the chosen coordinate system, the rotation from the z axis to the x-axis is positive. In 
Fig. 4, the W components of space velocity (the linear velocity along the z-coordinate) are 
plotted against heliocentric distance R. The points are plotted by using the data obtained 
by Drimmel et al. (2000) for a sample of 4250 Hipparcos OB stars in the magnitude interval 
0-13™. A comparison of Fig. 3 (the Sun is assumed to be at a Galactocentric distance of 
about Ro = 8.0 kpc) and Fig. 4 leads us to conclude that the value of = —0.36 mas yr“^ 
that we inferred from Hipparcos data almost completely explains the slope of the plot in 
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Fig. 4. Eliminating the dependence that we found should significantly reduce the precession 
of —25 km s“^ kpc“^ for distant OB stars derived by Drimmel et al. (2000). In this case, the 
fact that distant OB stars belong to the hydrogen layer implies that there is only a linear 
displacement of all distant OB stars along the z coordinate. 

CONCLUSIONS 

Based on a linear Ogorodnikov-Milne model, we have performed a kinematic analysis 
of the Hipparcos and TRC stellar proper motions. We found the rotation of all distant 
(r > 0.2 kpc) Hipparcos stars with a mean angular velocity of = —0.36 ± 0.09 mas 
yr“^ about the Galactic y axis. One of the causes of this rotation may be an inaccuracy 
of the lunisolar precession constant adopted when creating the ICRF (Ma et al. 1998). We 
showed that, in this case, the correction to the lAU (1976) constant of lunisolar precession 
in longitude is Api = —3.26 ± 0.10 mas yr“b We have shown that eliminating the rotation 
= —0.36 ± 0.09 mas yr“^ from the proper motions of distant OB stars should lead to 
a significant reduction in the -25 km s“^ kpc“^ precession of distant OB stars inferred by 
Drimmel et al. (2000). In this case, the fact that distant OB stars belong to the hydrogen layer 
reduces only to a linear displacement of all distant OB stars along the z coordinate. We have 
determined the mean Oort constants A = 14.9 ± 1.0 km s“^ kpc“^ and B = —10.8 ± 0.3 km 
s“^ kpc“^ from TRC data. The model component that describes the rotation of all TRC stars 
about the Galactic y axis was found to be nonzero for all magnitudes, = —0.86 ± 0.11 
mas yr“^. 

ACKNOWLEDGMENTS 

This work was supported by the Russian Foundation for Basic Research (proiect no. 
02-02-16570). 

REFERENCES 

1. T. A. Agekyan, B. A. Vorontsov-Vel’yaminov, V. G. Gorbatskii, et ah, A Course on 
Astrophysics and Stellar Astronomy, Ed. by A. A. Mikhailov (Fizmatlit, Moscow, 1962), Vol. 
2 [in Russian]. 

2. V. V. Bobylev, JOURNEES 1997, Ed. by J. Vondrak and N. Capitaine (Obs. de Paris, 
Paris, 1997), p. 91. 

3. V. V. Bobylev, Pis’ma Astron. Zh. 30, 185 (2004). 

4. B. W. Carroll and D. A. Ostly, An Introduction to Modern Astrophysics (Addison- 
Wesley, New York, 1996). 

5. P. Chariot, O. J. Sovers, J. G. Williams, et al., Astron. J. 109, 418 (1995). 

6. S. V. M. Clube, Mon. Not. R. Astron. Soc. 159, 289 (1972). 

7. S. V. M. Clube, Mon. Not. R. Astron. Soc. 161, 445 (1973). 

8. W. Dehnen and J. J. Binney, Mon. Not. R. Astron. Soc. 298, 387 (1998). 

9. R. Drimmel, R. L. Smart, and M. G. Lattanzi, Astron. Astrophys. 354, 67 (2000). 

10. E. H0g, A. Kuzmin, U. Bastian, et ah, Astron.Astrophys. 333, L65 (1998). 

11. A. Kuzmin, E. Hpg, U. Bastian, et ah, Astron.Astrophys., Suppl. Ser. 136, 491 (1999). 

12. P. G. Kulikovskii, Stellar Astronomy (Nauka, Moscow, 1985) [in Russian]. 

13. C. Ma, E. F. Arias, T. M. Eubanks, et ah, Astron. J. 116, 516 (1998). 

14. M. Miyamoto and M. Soma, Astron. J. 105, 691 (1993). 

15. B. du Mont, Astron. Astrophys. 61, 127 (1977). 


9 



16. B. du Mont, Astron. Astrophys. 66, 441 (1978). 

17. A. F. G. Moffat, S. V. Marchenko, W. Seggewiss, et al, Astron. Astrophys. 331, 949 
(1998). 

18. K. F. Ogorodnikov, Dynamics of Stellar Systems (Fizmatgiz, Moscow, 1965). 

19. R. P. Oiling and W. Dehnen, Bull. Am. Astron. Soc. 31, 1379 (1999). 

20. M. A.C. Perriman, L. Lindegren, J. Kovalevsky, et ah, Eur. Space Agency 1-17 
(1997). 

21. S. P. Rybka, Kinemat. Fiz. Neb. Tel 11, 77 (1995). 

22. R. L. Smart, R. Drimmel, M. G. Lattanzi, et ah, JOURNEES 1997, Ed. by J. Vondrak 
and N. Gapitaine (Obs. de Paris, Paris, 1997), p. 179. 

23. The HIPPARGOS and Tycho Catalogues, ESA SP-1200 (1997). 

24. V. V. Vityazev, Doctoral Dissertation in Mathematical Physics (St. Petersburg, 1999). 

25. H. G. Walter and C. Ma, Astron. Astrophys. 284, 1000 (1994). 


10 



